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ABSTRACT 

We study the stability of spectral approximations to scalar hyperbolic initial-boundary 
value problems with variable coefficients. Time is discretized by explicit multi-level or Runge- 
Kutta methods of order < 3 (forward Euler time differencing is included), and we study 
spatial discretizations by spectral and pseudospectral approximations associated with the 
general family of Jacobi polynomials. We prove that these fully explicit spectral approx- 
imations are stable provided their time-step, At, is restricted by the CFL-like condition, 
At < Const. AT -2 , where N equals the spatial number of degrees of freedom. We give two in- 
dependent proofs of this result, depending on two different choices of appropriate L -weighted 
norms. In both approaches, the proofs hinge on a certain inverse inequality interesting for 
its own sake. Our result confirms the commonly held belief that the above CFL stability 
restriction, which is extensively used in practical implementations, guarantees the stability 
(and hence the convergence) of fully-explicit spectral approximations in the non-periodic 
case. 
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1. INTRODUCTION. 


We are concerned here with fully-discrete spectral and pseudospectral approximations to 
scalar hyperbolic equations. In this context, the spectral (and respectively, the pseudospec- 
tral) approximations consist of truncation (and respectively, collocation) of N - term spatial 
expansions, which are expressed in terms of general Jacobi polynomials; Chebyshev and Leg- 
endre expansions are the ones most frequently found in practice. In this paper we prove that 
such iV-terms approximations are stable, provided their time-step, At, fulfills the CFL-like 
condition, A t < Const. N~ 2 . To clarify the origin of such CFL-like conditions in our case, 
we recall that the Jacobi polynomials are in fact the eigenfunctions of second-order singular 
Sturm-Liouville problems. Our arguments show that the main reason for the above CFL 
limitation is the 0(N 2 ) growth of the corresponding W-eigenvalues associated with these 
Sturm-Liouville problems. 

The paper is organized as follows. Section 2 includes a brief summary on the properties 
of Jacobi polynomials (and their quadrature rules) which are used throughout the paper. 
In Section 3 we state our main stability theorems for Forward Euler time-differencing and 
(pseudo-)spectral spatial differencing, for constant coefficients equations with homogeneous 
boundary conditions. Section 4 extends our stability results to the inhomogeneous case. In 
Section 5 we discuss multi-level and Runge-Kutta time differencing. Finally, in Section 6 we 
show how to extend our results in the presence of (positive) variable coefficients. 

2. VERY SHORT GUIDE TO JACOBI POLYNOMIALS. 

Jacobi polynomials, P* j3 , are the eigenfunctions of the singular Sturm-Liouville problem 
(2.1a) = -1 < a < 1, 

with corresponding eigenvalues A*, 

(2.16) A k = Afc(a, /?) = h(k + « + /? + 1). 

Different families of Jacobi polynomials are associated with different weight functions w(x), 
(2.1c) tu(x) = tu(x; a,/3) = (1 — x)“(l + *)^» a ,(3>— 1. 

In the sequel we shall frequently use several properties of the Jacobi polynomials. A brief 
summary of these properties is given below (consult e.g. [13]). We start with the well-known 

PROPERTY 1 ( Orthogonality) . We have 

(2.2) (Pf S ,Pi“ J ’ ) U,) = 0, i / *• 
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The derivatives of Jacobi polynomials are also Jacobi polynomials. This is evident from 
the following property which shows that {P* Q ’^ }*><> are orthogonal with respect to the 
weight (1 — x 2 )w(x) = w(x ; a + l,j3 + 1) and hence 

(2.3) = Const fc>a ^P^ Q+1,/9+1) , Constfc >a> /9 = ^(k + a + + 2). 

PROPERTY 2 ( Orthogonality of derivatives ). We have 

( 2 - 4 ) (Pj a ‘^ , P^' 1 = 0, j ^ h, 

( 25 ) 

Indeed, (2.4) and (2.5) follow from integration by parts of (2.1) against Pj a, ^\x). 


Let TTff denote the space of algebraic polynomials with degree < N. A useful consequence 
of the last two properties is provided by 

LEMMA 2.1 ( Inverse inequality) . For all per n we have 
( 2 - 6 ) ||p1I(i-x 2 Hx) < v/^IIpIUx), peiTff. 

Here u>(x) stands for an arbitrary w(x\ a, /3) weight, and \n = Ajv(a, (3) is the corresponding 
N-th eigenvalue. 


Remarks . 

1. The inequality (2.6) can be viewed as the algebraic analogue of the trigonometric 
inverse inequality, 

(2.7a) IIp'Hl’I-x,*] < (V||p|U 3 [-x.ir], P = any N - trigonometric polynomial. 

This should be contrasted with a similar L 2 -inverse inequality for algebraic polynomials 
where there is a loss of JV 2 -factor for each derivative [3], 

(2.76) l|p'IU 3 [-i,ij < Const. iV 2 ||p||z,j[_i i i], p = any N — algebraic polynomial, 

/_ i _ i \ 

and this estimate, (2.7b), is sharp in view of, e.g., p#(x) = J2k=o P* * K x )- Thus, the use 
of the different weighted L 2 -norms in the algebraic case, (2.6), is essential in order to retain 
a loss of only \/\n ~ N-iaciot for each derivative. 
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2. The inverse inequality (2.6) can be viewed as an i 2 -weighted version of Bernstein’s 
inequality 

||(1 - a 2 )V(*)IU~[-M] < JV|IK*)IU-[-i.il» P e7r "’ 

Standard interpolation arguments between this L°°-type estimate and the L 2 -type estimate 
(2.6) yield for q > 2 


(2.7c) ||(1 - * 2 )V(®)IU' ( . ;ai<J) M.i] ^ p€ ' KN ' 


Similar weighted L’-type estimates apply to higher derivatives. 


PROOF. Given p(x) in 7 r^, we will use its Jacobi expansion, p(x) — ^ A 1 ) 

and p'(x) = EfcLo a fc- P fc 0t, ^ ( x )- Starting with the left-hand side of (2.6) and using (2.4), (2.5) 

and (2.2) in this order, we obtain 


N 


(LHS) 2 = E«*H F ^ ), |l(i-« a M-) 


k= o 


= E < MRHsy. 


□ 


We note in passing that Lemma 2.1 can be generalized to higher derivatives: successive 
application of (2.6) with oj(x) = to(x; a, /?) yields 

(2.8) ||p ( * ) (®) ll(i— **)*«»{*) - II ^( a + + j) ' IIK z )llw(*)) P e *N- 

3 = 0 

This leads us to a ‘natural’ definition of non-periodic Sobolev spaces equipped with finite 
tf' (i) -norm, where, 

(2.9) IIpIIh^, = Sllp (fc) ll(i— *)*"(*)• 

With this in mind, we now recover a sharp inverse inequality familiar from the trigonometric 
setup 

(2.10) IIpII^j.) — Const, ‘ ^'||plk-)» Const ~ 1 + — , p€ir if. 

In the above discussion we can replace integrals by discrete summations in view of the 
well-known 

PROPERTY 3 ( Gauss quadrature rule) . Let {gAr(x)}^>o be a family of ir N -polynomials 
orthogonal w.r.t. the u(x)-weighted L 2 inner product. Let -1 < x x < x 2 . . . < x N < 1 be 
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the N zeroes of qn(x). Then there exists positive discrete weights , {wy}^!L lf such that for all 
peTT 2 N-i we have 

,1 N 

(2.11) / Lu(x)p(x)dx = ^2 Wjp^xf), pex 2 N-\- 


Remark . To compute the Gauss weights we set p(x) = j n (2.11). Since p(xf) = 0, j / k, 
(2.11) yields 

(2-12) u/ fc = — — - f u>(x) ?Ar ( Z ) dx, 1 < k < N. 

<lN( x k)J- 1 x ~X k ~ ~ 


PROOF. We have p(x) = t(x)qf/(x) + r(x) for some f(x) and r(x) in 7T N— \ . The choice 
of weights in (2.12) guarantees that (2.11) is valid for 7Tjv_i-polynomials, for 


span 


/tsfeU 

1 x - x k j 


= T if-1 • 


KJfcCJV 


This together with our assumption that qff^x} is -orthogonal to 7rjv" 1 imply 

yl i-l N N 

/ w(x)p(x)dx = / u>(x)r(x)<fx = J] ^j r { x j) = £ "iPC^i)- 

• / ~ 1 •'“l J--. 


J=1 


J=1 


EXAMPLES. 

1- Gauss-Jacobi quadrature rule . By Property 1, (2.11) applies to {-PaT^}a/>i with 
o;(x) = w(x;a,(3). Hence there exist {wj = wf(a,0)}^L 1 such that 

,1 N 

(2.13) / t u(x)p(x)dx = Wjp(xj ), for all pe7T 2 jv-i. 

i ._i 


Remark . The Gauss-Jacobi quadrature rule (2.13) can be used as a highly accurate quadra- 
ture rule for general smooth, not necessarily polynomial functions. The error incurred in 
such cases is governed by [4, p. 75] 

r 1 N 

(2.14) / w(x)f(x) -Yl w if( x i) = Const. / (2W) (0), Const. >0, |0| < 1. 

J ~ x i= l 
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2. Gauss-Lobatto-Jacobi quadrature rule . By Property 2, (2.11) applies to }w>o 

with oj(x) = (1 — x 2 )w(x; a, j3), and therefore, there exist {v)j = wf(a + 1 , /3 -f l)}£Li such 
that 


(2.15) 



N 

— x 2 )w{x)r(x)dx = ^2 ^i r ( x j)> f° r a U re7r 2 AT-l- 
3 = 1 


This is in fact a special case of the Gauss-Lobatto-Jacobi quadrature rule which is exact for 
all pe^ 2 N+i- Indeed, for all such p’s we have p(x) = (1 — x 2 )r(x) + £(x) with r(x) in TT 2 N -1 
and a linear l{x) = p(-l) 1 ^ + p(l)H £ - B V ( 2 - 15 ) 

Sli w(x)p(x)dx = Ejli Wjr(xj) + /i, w(x)t(x) = 

= T.U !%?<*,) + - x&. = i+u+ in- 


Thus, we have 


I = D w j K x i), a = 1 ^ - 2 ( a + P + *) 

4 ■*> I "'j 


J=1 


and the two expressions, II + III, amount to a linear combination of p(— 1) and p(l) 


II + III = WqP(x o) -I- 10 ^ +1 p(xjv + i), Xo = — 1 < Xi < . . . < Xjf < 1 = Xff+ 1 - 
Hence, there exist {wj = w^(a, such that 

f 1 

(2.16) / ti/(x)p(x)dx = 2^ iw^Xj), for all pe7r 2 j\r + i. □ 

- 7 - 1 j = o 


Finally, we shall need some information on the behavior of the collocation points which 
appear on the right of (2.13) and (2.16). We have 


PROPERTY 4 ( Distribution of zeros b If Xj = cosOj is the j-th zero of P^'^^x), then 
[1, p. 287] Ndj is a corresponding zero of Bessel’s function, and hence 


1 - x) = sin 2 ~ Const i0 W- 2 for jeJ = {!<;'< j 0 , N - j 0 < j < N }. 


Thus, the zeros of Pjy*’^(x) are accumulated within 0(N 2 ) neighborhood of { — 1, +1}- More 
precise estimates, e.g., [12, p. 19] yield 


(2.17) 


1 W 2 

1 — Xj ~~ 2(1 + a )’ 


1 <j <N. 
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3. FORWARD EULER WITH HOMOGENEOUS BOUNDARY CONDITIONS. 

We start with the scalar constant coefficient hyperbolic equation, 

(3.1) u t = au x , (x,t)e[— 1,1] x [0,oo), a > 0, 

which is augmented with the homogeneous condition at the inflow boundary, 

(3.2) tt(l, t) = 0, <>0. 

To approximate (3.1), we use the forward Euler time differencing on the left, and either 
spectral or pseudospectral differencing on the right. Thus, we seek a temporal sequence of 
spatial XAf-polynomials, v m = v N (x,t m = mAt), such that 

(3.3a) vjv(x, t m + At) = v N (x, t m ) -f At ■ v' N (x, t m ) + At • T(t m )q N (x). 

Here, qtf(x) is a 7r^r- polynomial which characterizes the specific (pseudo)spectral method we 
employ, and r = r(t m ) is a free scalar multiplier to be determined by the boundary constraint 

(3.36) v N (x = 1, t m ) = 0. 

We shall study the spectral- Jacobi tau methods, [8], [2], where 

(3.4) v N (x, t m + At) = v N (x, t m ) + At • av' N (x, t m ) + At • T(t m )gjv(x), qti{x) = Pfc' 0 \x), 
and the pseudospectral-Jacobi methods, [5], [2], which are collocated at the interior extrema 

n f p( a <0) : p 

(3.5) ujv(x, t m + At) = Vff(x, t m ) + At • av'pf(x, t m ) + At • r(t m )^(x), ?at(x) = Pj?+i (*)• 

Remark. These two families of spectral and pseudospectral Jacobi methods are closely related 
since (x) is a scalar multiple of Pjy* +1 ’^ +1 ^(x), consult (2.3). We will not discuss here a 

different alternative to (3.5) where one collocates at the interior extrema of P^' 0 \x) together 
with the downstream outflow boundary so that 

(3.6) qw(x) = (1 + x)P { “' 0) (x). 

Let —I < Xi < x 2 < . . . < xn < 1 be the ^-different zeros of the forcing polynomial 
The spectral approximation (3.3a) restricted to these points reads 

(3.7a) v N (xj, f m+1 ) = v N {xj, t m ) + At • av' N (xj, t m ), 1 < j < N, 
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and is augmented with the homogeneous boundary condition 


(3.76) t m ) — 0. 

Equations (3.7a), (3.7b) furnish a complete equivalent formulation of the spectral approxi- 
mation (3.3a), (3.3b). An essential ingredient in a stability theory of such approximations 
lies in the choice of appropriate L 2 - weighted norms 

N 

(3.8) ll/MIIL =</(*)./(*)>. 

j=l 


We make 


DEFINITION 3.1 ( Stability) . The approximation (3.7a), (3.7b) is stable if there exist 
discrete weights, {uj, > 0}^, and a constant rj 0 independent of N> such that 

(3.9) ||vn(*>*)||« ^ Const. e’’ ot ||v JV (-,0)||„. 

The approximation (3.7a), (3.7b) is strongly stable if (3.9) holds with Const. = 1 and rjo ^ 0, 
i.e., if 

(3.10) ||vtf(',t)|| u < ||vAr(-,0)|| u ,. 

We recall that in the Jacobi-type spectral approximations (3.4) and (3.5), the nodes 
are the zeros of Jacobi polynomials associated with the Gauss and Gauss-Lobatto 
quadrature rules. We use 

(3.11) — min(l + Xi, 1 — Xff) 

to measure the minimal gridsize associated with these Gauss nodes. Our choice of discrete 
weights {uj}jLi for the stability of the spectral and pseudospectral- Jacobi methods (3.4), 
(3.5) will be specified later on; these weights are related (but not equal) to the corresponding 
Gauss weights { Wj }^ =1 indicated earlier. 

With this in mind we have 


THEOREM 3.1 ( Stability of the spectral and pseudospectral Jacobi methods) . Con- 
sider the spectral approximations (3.7a), (3.7b), associated with the Jacobi-tau method (3.4) 
or the pseudospectral- Jacobi method (3.5). There exists a positive constant T) 0 = Vo(&, P) > 0 
independent of N, such that if the following CFL condition holds 


(3.12) 


At • a 



—) 


< Vo, 
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then the approximation (3.7a), (3.7b) is strongly stable and the following estimate is fulfilled 


(3.13) 




Notes. 


1. The choice of L 2 -weighted norms . Theorem 3.1 deals with the stability of both the 
spectral-tau methods associated with Jacobi polynomials a,/?e(— 1,0), and with 

the closely related pseudospectral methods associated with ■P/v+i^ ( x )> a:, — 1,0). In each 

case we give two different stability proofs which are based on two different choices of discrete 
Zr 2 -weighted norms; these discrete weights are given by 


(3.14a) 




1 + gj 

1 — Xj 


Wj, where = Gauss — Jacobi weights in (2.13), 


(3.14ft) u>j = (1 + Xj)wj } where { Wj }yL x = Gauss — Lobatto Jacobi weights in (2.16). 


2. The CFL condition . The CFL condition (3.12) places an 0(W~ 2 ) stability restriction 
on the time-step At, and this stability restriction involves two factors. First, since we expand 
our solution in terms of the eigenfunctions of the Sturm-Liouville problem (2.1), the CFL 
condition involves the corresponding ( N — l)-th eigenvalue 


(3.15a) 


= A^_i(a, /3) < N 2 , a,/3e(— 1,1). 


Second, the spectral Jacobi approximation (3.7a) is collocated at Gauss nodes which are 
accummulated within 0(N ~ 2 ) neighborhoods near the boundaries, i.e., by (2.17), 


(3.156) 


1 N 2 

AXmin “ 2(1 + a)' 


Thus, in view of (3.15a), (3.15b), the CFL condition (3.12) boils down to 


(3.16) 


At • aN 2 < Const., 


Const. = 


Vo 


1 + a 

2-fa 


> 0. 


3. The choice of a stability norm . The stability statement asserted in Theorem 3.1 is 
formulated in terms of discrete semi-norms, || • H^, which are u - weighted by either (3.14a) or 
(3.14b). We note that || • are in fact well-defined norms on the space of 7T^-polynomials 
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satisfying the vanishing boundary condition (3.7b), i.e., corresponding to (3.14a) or (3.14b) 
we have 

(3.17a) |M-,*)IU = j\ w(x)^^v 2 N (x } t)dx, v N (l ,i) = 0, 

(3.176) \\v N {-, OIL = + s ) v n( x > *)<**> "M 1 * 0 = °* 


Moreover, 
usual L 2 -type 


in view of (3.15b), one may convert the stability statement (3.13) into the 
stability estimate at the expense of possible algebraic growth which reads 


(3.18) |M-,t)||«(«) < ,wat !M->0)IL ( x ) , IM'.OII^x) - J_^w(x)v 2 N (x,t)di 


4. Exponential time decay . Let us integrate by parts the differential equation (3.1) 
against (1 + x)u. Thanks to the homogeneous boundary condition (3.2) we find 

( 3 . 19 ) ^(1 + x)u 2 (x,t)dx < -~J{\-\-x)u 2 {x,t)dx, 

and therefore, 

(3.20) IW*i0lli+* < e_ * ot || u ('» 0)||i+** 

This estimate corresponds to the special case of the stability statement (3.13) for the spectral- 
Legendre tau method (« = /? = 0) weighted by (3.14b). The exponential time decay indi- 
cated in (3.20) and more generally in (3.13), is due to the special choice of unweighted 
stability norms. The weights {wj}^ in (3.14a), (3.14b) involve the essential factors 1 + x, 
or l±|i which amplify the inflow boundary values in comparison to the outflow ones. Since in 
the current homogeneous case, vanishing inflow data is propagating into the domain, this re- 
sults in the exponential time decay indicated in (3.20) and likewise in the stability statement 

(3.13) . 

5. The inflow problem. A stability statement similar to Theorem 3.1 is valid in the inflow 
case where a < 0. Assume that the CFL condition (3.12) holds with r, 0 = Vo{P,<*), then 

(3.13) follows with discrete weights Uj = or = C 1 “ Vj) w j- 

The rest of this section is devoted to the proof of Theorem 3.1 according to the various 
cases outlined in the four lemmata below. We start with 
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LEMMA 3.2 ( Stability of the spectral-tau method) . Consider the spectral- Jacobi-tau 
method (3.4). Then Theorem 3.1 holds with 

(3.21a) uij = - — ~ w j, {u>j = w f(<x,P)}jLi = Gauss — Jacobi weights, 


(3.216) 


Vo = Vo (<*,/?) 


' a+/3< 0 , 

. |(1 -«), a + /3> 0. 


a,/?e(-l,l). 


PROOF. Squaring of (3.7a) yields 

IM-,* m+1 )l£ = IM-,r)l£+ 

( 3 - 22 ) +2 At ■ a < v N (-,t m ),v' N (-,t m ) > +(At • a) 2 ||u^(-,t"*)|| 2 = 

= IM-, t m )\\l + 2A t • al + (A t ■ a) 2 1 1 , 

and we turn to estimate the two expressions, I and II, on the right of (3.22). 

First let us note that since the 7r //-polynomial f//(x, t m ) vanishes at the inflow boundary, 
(3.3b), we have 


( 3 - 23 ) v N (x, < m ) = (1 - x)p{ x) for some p(x ) = p^-iO^ex/v-i. 


Also, a straightforward computation shows that 

(3.24) ^(x)——^ (i _ x f = [(/3 - a + 2) - (/? + a)x]u;(x) > 4t] 0 w(x), \x\ < 1, 

where tj 0 = rj 0 (a,f3 ) is given in (3.21b). 

Now, since f^fujv(x, i m )v' N (x, t m )eTr 2 N-i> Gauss quadrature rule (2.11) implies 

I = w j7 _ J v N(xj, t m )v' N (xj, t m ) = f -uj(x)_t — v N (x, t m )v' N (x, t m )dx. 
j = i 1 *'—1 1 — x 

We integrate by parts the RHS of I, substitute v N (x, t m ) = (1 - x)p(x) from (3.23), and in 
view of (3.24) we o btain 

(3.25) / = — i /* (™(x)i±^)'(l - x)V(x)*: < -2 % || p ||^,. 


Next, let us consider the second expression, II, on the right of (3.22). As before, we substitute 
vjv(x,f n ) = (1 — x)p(x) from (3.23) and obtain 


N 


II = mil = - Xi)p'(xy) - P(x,)] ! < 


i = 1 


U 1 

3 1 -x 
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N 


1 + Xj 


< 2 £ wjl - x?)(p'(*;)) 2 + 2 E "ir-T-P M = 11 1 + 7/ ” 

;=i i=i 1 

Since (1 - i 2 )(p # (i)) 2 e 7 r 2 ^- 2 ) the Gauss quadrature rule (2.13) followed by the inverse in- 
equality (2.6) implies 

II\ = 2^w i (l - ®j)(p'( s j')) 2 = 2 IIp , |I(i-* j Hx) ^ 2Ajv-i||p||L(*)> P^jv- 

i= i 

This together with the obvious upper-bound 

* 1 + 


Hi = 


< Tr-llpllt(x), 


gives us 
(3.26) 


II 


— ^2Ajv_i + IIpIIw(x)- 


Equipped with (3.25) and (3.26) we return to (3.22) to find 

2t/ 0 — A t • a ^Ajy_i 


+ 


Ax r 


-)] IWIlw- 


(3.27) \\v H {; t m+1 ) \\l < |M*. Oil" “ 2Ai • a 

The CFL condition (3.21b) implies that the squared brackets on the right are nonnegative, 

AZr)]^” >0 ' 


(3.28) 


2t}q — A t 


o. ^Aw_i 


+ 


and hence strong stability. 

In fact, one more application of the Gauss quadrature rule yields 


(3.29) 




mMI2 


4 mM 2 
u ; 


The inequalities (3.29), (3.28) together with (3.27) imply 
(3.30) |K(-,t m+1 )||?, < (1 - 2r7 0 At- a)||vtf(-,t m )|| 

and the result (3.13) follows. ^ 

Next, we take advantage of the rather general setup we employed in Lemma 3.2. Specif- 
ically, since is proportional to Pfr +l ’ 0+1) , consult (2.3), we may use Lemma 3.2 with 

Tjo(a,(3 ) replaced by rj 0 (a + 1,/? + 1) to conclude: 
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LEMMA 3.3 ( Stability of the pssudospectral method ). Consider the pseudospectral- 
Jacobi method (3.5). Then Theorem 3.1 holds with 

(3.31a) ujj = — - -Wj y + l,/3 + 1)1 = Gauss — Jacobi weights, 

1 ' j — 1 


(3.31.6) ri Q = r) 0 (a,/3) = -^>0, a,/3e(- 1,0). 

As mentioned before, alternative proofs of Theorem 3.1 are possible. For example, fol- 
lowing [6, Theorem 5.1] one may employ a stable norm weighted by uij = (1 + Xj)wj instead 
of the u)j = jrfHoj weights used in (3.21a), (3.31a). We begin with 


LEMMA 3.4 ( Stability of the spectral-tau method revisited) . Consider the spectral- 
Jacobi tau method (3.4). Then Theorem 3.1 holds with 

(3.32a) u>j = (1 + Xj)wj, { Wj = u»®(a, )S)J = Gauss — Jacobi weights 


(3.326) 


Vo = Vo (a, /3) 


< 


a 
2 » 


i K 1 - 


Ct + p + 1 > o, 

a,/3e(— 1,0). 

a + /? + 1 < 0, 


PROOF. We square (3.7a) - this time using the inner product weighted by (3.32a), 

|M-,t m+1 )|£ = IM-.OII* +2At.a<UAr(- > t m ),<(.,t") > + 

(3.33) +(At.a) 2 |K(-,mii = 

= ||»/»r(-7^)|li + 2 At • a/ + (At • a) 2 II, 

and as before we have to estimate the two expressions on the right of (3.33) 

The first expression, I, involves discrete summation of f(x ) = (1 + x)v N {x, t m )v' N {x, f m ), 
and since /(x) is a 7T2iV"Polynonual, its N - nodes Gaussian sum is not an exact integral. The 
essential observation here is that / (2A ^ = Const. > 0 in this case, and the error estimate 
(2.14) tells us that the Gauss quadrature rule is upper bounded by 

N ,i 

(3.34) 7 = 2 Wj{\ + Xj)v N (xj, t m )v' N (xj, t m ) < / ™(x)(l + x)v N (x, t n )v' N (x, t m )dx. 

3 = 1 J ~ 1 

Let us recall that the homogeneous inflow boundary condition (3.7b) implies 

(3.35) VN(x,t m ) = (1 — x)p(x), for some p(x) = pn-i{x)e'KN-i. 
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Also, a straightforward computation shows 

(3.36) (u>(x)(l + x))'(l -x) = [(/?-a + l)-(a + /3 + l)x]u>(x) > irj 0 w{x), |x| < 1 
where r} 0 = 7/ 0 (a,/?) is given in (3.32b). 

We integrate by parts the RHS of I, substitute (3.35), and in view of (3.36) we obtain 

(3.37) I = — 1 J {w{x){ 1 + x))'(l - x) 2 p 2 (x)dx < -27/ 0 ||p||(i_ I ) u ,( I )- 

Next, let us consider the second expression, II, on the right of (3.33). As before, we substitute 
Vff(x,t m ) = (1 - x)p(x) from (3.35) and Gauss quadrature yields 

II = E£i «i(l + *,-)[(! - *M*i) - K^-)] 2 = 

(3.38) = Sli «>(*)(1 - ® 2 )(1 “ x)(p'(x)) 2 - 2 /-i w (*)( 1 “ x 2 )p(x)p’(x)dx 

+ u;(x)(l + x)p 2 (x)dx = IIi + H% + Us- 

The inverse inequality (2.6) with weight u>(x) = (1 — x)u?(x) implies 

Hi = ||p , ||(i-**X 1 -*M*) — A ^-iIIpII(i-xMx)» ^n-i = ^N- i( a + l>/5)> 

and integration by parts of II2 gives 

Ih + Ih = /-i[K*)(l - * 2 ))' + w(x)(l + x)]p 2 (x)dx 

< 2||p 2 |U (x) = 2Y.U w^/ixi) < ^-||p||?x_ sMx) . 

Consequently we have 

(3.39) II ^ ^AT— 1 “I" “j ||p!I(1— x)w(x)' 

Equipped with (3.37) and (3.39) we return to (3.33) to find 

(3.40) |M-,< m+1 )ll« < IM-.OII 2 “ 2A t a [2* - At • a (a w _ x + ^-)] IIpII(i-*M*) 

and the result follows along the lines of Lemma 3.2, consult (3.27). ^ 

Lemma 3.4 does not cover the pseudospectral Jacobi methods, since by (2.6) the corre- 
sponding Jacobi parameters a + l,/3 + 1 ^ (—1,0). However, the proof of Lemma 3.4 can 
be carried out in the pseudospectral case if we replace the Gauss quadrature rule by the 
Gauss-Lobatto one. We omit the almost identical details (which are outlined for the variable 
coefficients case in Theorem 6.2 below) and state 
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LEMMA 3.5 ( Stability of the pseudospectral method revisited) . Consider the pseu- 
dospectral Jacobi method (3.5). Then Theorem 3.1 holds with 

(3.41a) ufj = (1 + Xj)wj, {wj = Wj(a, l3)}y =1 = Gauss — Lobatto — Jacobi weights, 

f ~b a + f3 + 1 > 0, 

(3.415) rjo = T) 0 (a, (3) = l a,/?e(— 1,0). 

[ a + /3+l <0. 

Remark. The stability asserted in Lemma 3.5 is stated in terms of the discrete semi- 
norm ||v^f(-, t)||* = w j(l + x j) v ff( x j, t) weighted by the interior Gauss-Lobatto weights 
However, taking into account the homogeneous boundary condition (3.7b) and the 
exactness of Gauss-Lobatto quadrature for tt 2 N+i polynomials, this amounts to 

N+l 

II v aK-.*)II« = £ w j { 1 + x j) v N( x i, *) = 

J=0 

4. FORWARD EULER WITH INHOMOGENEOUS INITIAL-BOUNDARY 
CONDITIONS. 

We consider the inhomogeneous scalar hyperbolic equation 
(4-1) u t = au x + F(x, t), (x,f)e[— 1, 1] x [0, oo), a > 0, 

which is augmented with inhomogeneous data prescribed at the inflow boundary 
(4.2) <1 ,t) = g(t), t> 0. 

Using the forward Euler time differencing, the spectral approximation of (4.1) reads, at the 
iV-interior zeros of q^(x) } 

(4-3 a) v N (x j ,f n+1 ) = v N { Xj , t m ) + At ■ av' N {x h t m ) + A tF( Xj , t m ), q N ( Xj ) = 0, 

and is augmented with the boundary condition 

(4.34) Ml 

In this section, we study the stability of (4.3), (4.4) in the two cases of 
(4.4a) Spectral — Jacobi tau method : ^(x) = Pjf’^\x), a, (3e( — 1,0), 


J tt>(x)(l + x)vjy(x, t)dx. 
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and the closely related 

(4.46) pseudospectral — Jacobi method : Qn{x) = ( x )i a » P e i ~ ^)* 

To deal with the inhomogeneity of the boundary condition (4.3b), we employ a device 
introduced in [6, Section 5]. Namely we consider the -k n - polynomial 

(4.5) v N (x, t ) = v N (x, t) - ^yjsC*)- 
If we set 

(4.6) F(x,t) = F(x,t) + a^g(t) 
then V N (x,t ) satisfies the inhomogeneous equation 

(4.7a) V N (x j} t m+1 ) = V N (xj, t m ) + At ■ aV^Xj, t m ) + AtF( Xj , t m ), 

which is now augmented by the homogeneous boundary condition 


( 4 . 76 ) Vjv(M m ) = 0 . 

Theorem 3.1 together with DuhammePs principle provides us with a prion estimate of 
||Vftr(-,t)|U in terms of the initial and the inhomogeneous data, \\V N (-, 0)|| w and ||F(-,i)||». 
Namely, if the CFL condition (3.12) holds, then we have 

( 4 . 8 ) IIWC-.OIU ^ e“*°‘l|Vtf(->°)llu' + X) At-e ”° o(t * ) ||^’(-,* )IU 

o<t m <t 

Since the discrete norm || • ||„ is supported at the interior zeros of qff(x) where V N ( X j, t) = 
t) we conclude 


THEOREM 4.1 ( Stability of the spectral tau and pseudospectral Jac obi methods). Con- 
sider the spectral approximation (3.3a), (3.3b) associated with the Jacobi-tau method (4.4a) or 
the pseudospectral- Jacobi method (4.4b). There exists a positive constant rj 0 = r}o(a,P) > 0 
independent ofN, such that if the following CFL condition holds (consult (3.12)), 

( 4 . 9 ) Ai • a (a„_, + ) < Vo 


then the approximation (3.3a), (3.3b) satisfies the stability estimate 


(4.10) 

IM-,t)||„ < 


M-. 0 )H.+ E 




IkW- 


i?w(i)i 


%(OI 
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Theorem 4.1 provides us with a priori stability estimate in terms of the initial data, 
the inhomogeneous data, F(-,t), and the boundary data ^(<). The dependence on 
the boundary data involves the factor of which grows linearily with N, as shown by 

LEMMA 4.2. Let be the discrete weights given by either (3.32a) or (3.41a). 

Then there exists a constant independent of N such that 

( 411 ) < Const. N. 

Mi)| 


Remark . The stability estimate (4.10) together with (4.11) imply 
(4.12) 

ll^(-,<)|U < e -T,oat ||uj V (-,0)|| u ,+ At-e -TK)0 <‘- tm ) [||F(-, t m )|j w + Const.WAt • e - ’ ,0 ^ t - <m )| 5 (t m )|] 

An inequality similar to (4.12) is encountered in the stability study of finite-difference ap- 
proximations to mixed initial-boundary hyperbolic systems [9], We note in passing that the 
stability estimate (4.12) together with the usual consistency requirement guarantee the spec- 
trally accurate convergence of the spectral approximation; consult [7] for the semi-discrete 
case. : 


(a^) R °° F ‘ WC COnsider ’ for exam P le » the spectral-tau method associated with q N (x) = 
P" * and with dis ^ete weights = (1 + Xj )wf (a,/3). Using Gauss and Gauss-Lobatto 
quadrature rules (2.13) and (2.16) in this order, we obtain 

M-OC = £ «<*X1 + xXpfrn'MYdx = 2 »‘ +1 (a./s)|pt“’«'(i)|>, 

and (4.11) follows in view of 


2 w n+i(oi, 0) 


pfr®' to 




< Const. N 2 . 


Similar arguments (which are omitted) apply to the pseudospectral approximation associated 
with q N (x ) = (i) and with discrete weights u jj = (1 + Xj)wf(a,/3). 

5. MULTI-LEVEL AND RUNGE-KUTTA TIME DIFFERENCING. 


In the previous sections we proved the stability of spectral approximations which are 
combined with the first order accurate forward Euler time differencing. In this section we 
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extend our stability result for certain second and third order accurate multi-level and Runge- 
Kutta time differencing, which were studied in [10], [11]. 

To this end we view our 7 r^-approximate solution at time level t, v(-,t), as an (N + 
l)-dimensional column vector which is uniquely realized at the Gauss collocation nodes 

(v(xi, 

v(x N ,t),v(l,t)). 

The forward Euler time differencing with homogeneous boundary conditions (3.4), (3.5), 
reads 

(5.1a) v(t m + At) = [I + At • aL]v(t m ), a > 0, 

where L is an (N + 1) x ( N -f 1) matrix which accounts for the spatial spectral differencing 
together with the homogeneous boundary conditions 

(5.16) Lv(t m ) = (u'(xi, t m ), . . . , v'(x N , f m ), 0). 

Remark . The construction of a spectral differentiation matrix L can be accomplished m one 
of two ways. One possibilty is carried out in the physical space, by exact differentiation 
of the unique Tr^-interpolant which assumes the gridvalues v(x 1 , t), ..., v(x N , t), u(l, t) at the 
corresponding Gauss nodes. This leads to full (N + 1) x (N + 1) differentiation matrices 
L, which are recorded for example in [2]. Spectral differencing in the physical space is 
then carried out by a matrix-vector multiplication at the expense of 0(N 2 ) operations. 
Alternatively, spectral differentiation can be accomplished in the transformed space. In the 
particular case of Chebyshev (pseudo-)spectral method, this leads to a factorization of the 
corresponding differentiation matrix L, whose matrix-vector multiplication can be carried 
out efficiently by FFT requiring O(NlogN) operations, consult [8], [2]. 

Theorem 3.1 tells us that if the CFL condition (3.12) holds, i.e., if 

(5.2) A< ■ « (A«-i + i Vo, 

then I + At ■ clL is bounded in the cj-weighted induced operator norm, 

(5.3) \\I + At-aL\\ u <e-^ a - At . 

Let us consider an (s + 2)-level time differencing method of the form 

a 

(5.4) v(t m + A t)=J2 e k [I + CfcAf • aL]v(t m - k ), c k >0, 6 k > 0, ^ 0 k = 1. 

fc =0 k ~° 

In this case, v(t m + At) is given by a convex combination of stable forward Euler differencing 
and we conclude 
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COROLLARY 5.1 ( Multi-level time differencing) . Assume that the following CFL con- 
dition holds, 

(5.5) At • a 1 ^ Cfc > 0, k = 0, 1, . . . , s. 

\ LAX min' Cjfe 

Then the spectral approximation (5.4) is strongly stable and the following estimate holds 

(5-6) ll' u Jv(-,t)||a; < e -r> * ot j|uAr(-, 0 )|| & ,, 77 . = min — > 0 . 

* Ck 

In [10], second and third order accurate multi-level time differencing methods of the 
positive type (5.4) were constructed. They take the particularly simple form 

(5.7) v(t m + At) = 6[I + c 0 At ■ aL]v(t m ) + (1 - Q)[I + c.At • aL]v{t m -), 

with positive coefficients given in Table 1. 



Second-order time differencing | 0 | cq 


4- level method (s = 2) 

5- level method (s = 3) 


Third-order time differencing 


5- level method (s = 3) 

6- level method (s = 4) 

7- level method (s = 5) 



Table 1. Multi-level methods. 

Similar arguments apply for Runge-Kutta time differencing methods. In this case the 
resulting positive type Runge-Kutta methods take the form 

(5.8a) u ( 1 )(r +1 ) = [/ + At • aL}v{t m \ 

(5.86) u«(r +1 ) = 0 k v(t m ) + (1 - e k )[I + AtaLjv^t^ 1 ), k = 2, ■■■,*, 


v(t m+1 ) = V (‘\t m+1 ). 
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We state 


COROLLARY 5.2 (Runge-Kutta time differencing). Assume that the CFL condition 
(3.12) holds. Then the spectral approximation (5.8a-c) with 0 < Ok < 1 strongly stable and 
the stability estimate (3.13) holds. 

In Table 2 we quote the preferred second and third-order choices of [11]. 


Second order time differencing 

0 2 

^3 

Two-step modified Euler (3 = 2) 

1 

2 

- 

Third order time differencing 



Three-step method (s = 3) 

3 

4 

1 

3 


Table 2. Runge-Kutta methods 


Remarks . 

1. The above results can be extended to include nonhomogeneous data as well. We omit 
the details. 

2. In [10], [11], higher order accurate (> 3) multi-level and Runge-Kutta time differ- 
encing schemes were constructed. In the present context (of constant coefficient spectral 
approximations), they amount to convex combinations of the stable forward Euler differ- 
encing I + At ■ aL and its adjoint I - At ■ aL. The above argument does not cover these 
cases, however, since in our case the stability of I ± At • aL is measured by different weighted 

norms. 


6. SCALAR EQUATIONS WITH VARIABLE COEFFICIENTS. 


We begin with 

EPILOGUE. When dealing with finite-difference approximations which are locally supported, 
i.e., finite difference schemes whose stencil occupies a finit e number of neighboring grid cells 
each of which is of size Ax, then one encounters the hyperbolic CFL stability restriction 


( 6 . 1 ) 


— | a | < Const. 
Ax 
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With this in mind, it is tempting to provide a heuristic justification for the stability of 
spectral methods, by arguing that a CFL stability restriction similar to (6.1) should hold. 
Namely, when Ax is replaced by the minimal gridsize, Ax mtn = min, |x J+1 - Xj\ - 0(N~ 2 ), 
then (6.1) leads us to 

(6.2) At - |ajjV 2 < Const. 

Although the final conclusion is correct (consult (3.16)), it is important to realize that 
this hand waving argument is not well-founded in the case of spectral methods. Indeed, 
since the spectral stencil occupies the whole interval (-1,1), spectral methods do not lend 
themselves to the stability analysis of locally supported finite-difference approximations. Of 
course, by the same token this explains the existence of unconditionally stable fully implicit 
(and hence globally supported) finite difference approximations. 

As noted earlier, our stability proof (in Theorem 3.1) shows that the CFL condition (6.2) 
is related to the following two points: 

#1. The size of the corresponding Sturm-Liouville eigenvalues, A^-i = 0(N 2 ). 

#2. The minimal gridsize, ^ = 0(N 2 ). 

The second point seems to support the fact that Ax^ plays an essential role in the CFL 
stability restriction for the global spectral methods, as predicted by the local heuristic argu- 
ment outlined above. To clarify this issue we study in this section the stability of spectral 
a PP rox i ma ti°Jis to scalar hyperbolic equations with variable coefficients. The principle raison 
d etre which motivates our present study, is to show that our stability analysis in the con- 
stant coefficients case is versatile enough to deal with certain variable coefficients problems. 
We begin with the particular example introduced in [8] 

(6*3) u t — -xu X} (x, t)e[ — 1 , 1] x [0,oo). 

We shall show that the CFL stability restriction in this case is related to the 0(N 2 )- size 
of the Sturm-Liouville eigenvalues (point #1 above), but otherwise it is independent of the 
minimal gridsize mentioned in point #2 above. 

Observe that (6.3) requires no augmenting boundary conditions, since both boundaries, 
x = are outflow ones. Consequently, the various Trjv-spectral approximations of (6.3) 
with forward Euler time differencing, read 

(6-4) vn{x> + At) = v N (x , t m ) — At • xv' N (x y t m ). 

We have 
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THEOREM 6.1 ( Stability) . Assume that the following CFL condition holds 


(6.5) At-Aj\r<l, A# = N(N + 1). 

Then the spectral approximation (6.4) is strongly stable and the following estimate is fulfilled 

(6.6) [|vat(-, t)||i_*2 < ||ujv(-,0)||i_«j. 

PROOF. Squaring (6.4) yields 

(6.7) -2A<(v*(-,i m ),ZT/M-, O)!-*’ + (At) 2 ||xV^(-, t m )||i_ x a 

= \\v N (;t m )\\l_ x7 +2At-I + (At) 2 -IL 

Integration by parts shows that the first expression, I, is given by 

(6.8) I= l - J\x(l - x 2 ))'^(x, t m )dx = ||M-, t m )||?_ x i - J ^ x 2 v 2 N (x, t m )dx. 

Next, we write xv^f = (xu A r ) / — ujv, &nd by the inverse inequality (2.6) the second expression, 
II, does not exceed 

(6.9) II = ||(xv* (z, i m )Vi- x ’ ~ 2 £ * 2 «Sr(*» tm ) dx ^ ~ 2 ) £ x2v n( x > tm ) dx - 

Inserting (6.8) and (6.9) into (6.7) we end up with 

(6.10) |K(-,t m+1 )||?_ xJ < (1 + At) • ||v^(-, t m )||J_ I j + At- [(\ff -2)At — 2]- J ^ x 3 Vpf(x, t)dx. 
The CFL condition (6.5) tells us that the contribution of the second term is negative, for 

At[(Ajv - 2) At - 2] • £ x 2 v 2 N (x,t m )dx < -At||u(-, f m )||?_ x ’ 

and the strong stability estimate (6.6) now follows. D 

lAfe now turn to disscuss scalar hyperbolic equations with positive variable coefficients 

(6.11) u t = a(z)u x , 0 < a(x) < a*,, (z, t)e[-l, 1] x [0, oo), 

which are augmented with homogeneous conditions at the inflow boundary 

(6.12) u(l,t) = 0. 
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We consider the pseudospectral- Jacobi method collocated at the ./V-interior zeros of ( x ). 

Using forward Euler time differencing, the resulting approximation reads 


(6.13a) v N ( Xj , t m+1 ) = v N (xj, t m ) + At ■ a(x i K(z i> t m ), = 0, 

together with the boundary condition 

(6.136) ujv( 1, < m ) = 0. 

Arguing along the lines of Theorem 3.1, we have 


THEOREM 6.2 ( Stability of the pseudospectral Jacobi method with variable coefficients) . 
Consider the ■pseudo spectral- Jacobi approximation (6.13a), (6.13b). There exists a constant 
Vo = Vo{ a t P), 


(6.14a) 


Vo = Vo{a,P ) = < 


a 

' 2 ’ 


a+P+ 1 > 0, 


«,M-1,0), 


l |(1 - P), a + p + 1<0, 

such that if the following CFL condition holds 

(6.146) At ( a^N-i + 2 max j < Vo> 

then the approximation (6.13a), (6.13b) is strongly stable, i.e., there exist discrete weights 

\U ' f * 1 N 

(6.15a) u>j = ( 1 + Xj) . - . ■ , = w^(a, P)\ = Gauss — Lobatto weights, 


such that 


(6.156) 


IW',t)||ar < ||v w (-,0)||„. 


PROOF. We divide (6.13a) by <Ja(xj), 

- I = v N (xj ,t m+1 ) = -=L=v N (x jf t m ) + At ■ ^a(xj) • v' N (x, t m ), 

V“(*i) \l a { x i) 

and, proceeding as before, we square both sides to obtain 

IM-,<’" +1 )lll = IM-,f m )IIS+ 

(6.16) +2A t < w(-, 4”), >/„(•, t”) > +(A() J |M K(.,<")||J, 
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The first expression, I, involves discrete summation of the tc 2 n- polynomial f(x) = (1 + 
x)v N (x,t m )v' N (x,i m ) and since /(±1) = 0 (in view of (6.13b)), the ( N + l)-nodes Gauss- 
Lobatto quadrature rule yields 

N+l fl 

I = 53 w i t 1 + x i) v *( x v <”>*(*;> O = J + *)^(x, t m )v N ( x , t m )dx. 

3 = 0 

We integrate by parts the RHS of I, substitute v^(x,t m ) = (1 — x)p(x) from (3.35), and in 
view of (3.36) we obtain as before (compare (3.37)), 

(6.17) I < -2t7o||p|| (!_*)«,(*)• 


The second expression, II, gives us 

II = 1 ‘U'ja(zj)(l + Zj)[(l ~ x i)p'( x j) ~ K z i)] 2 ^ 

(6.18) < 2 aoo EjLx ™j(l + x i) 2 (p'( x j)) 2 + 2 £?* w ja (xj ) (1 + Xj)p 2 (xj) 

= 2aooII\ + 2 ■ Hi- 

The inverse inequality (2.6) with weight u;(x) = (1 — x)uj(x) implies 

Hi = IIp'II(1-* 3 )(1-*Mx) < *JV-l||p||(l-,) w (x)» = *n-i(<* + 1 >P) 


and the expression II 2 does not exceed 


III < 5 Wj ^ ~ x ^ p2 ( Xj } - 2 ' 1 - X } ' 

Consequently we have 

(6.19) II <2 ^ooAjv.x + 2 • max IIpII(i-xM*)- 


tf+1 


a ( x i) 


Equipped with (6.17) and (6.18) we return to (6.16) to find 

( 6 . 20 ) 


IM-,t m+1 )ll 2 <IM->OII 2 -2Ai 


2t?o 


— A t ^a 0 


«(*,-) 

+ 2 max 

l<j<Nl-Xj t 


IIpII( 


1 — x)ti/(x) J 


and (6.18b) follows in view of the CFL condition (6.14b). 


□ 


Notes . 

1. The case a(xj) = a = Const. > 0, corresponds to the stability statement of Lemma 
3.5. Similar stability statements with the appropriate weights which correspond to Lemmata 
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3.2, 3.3, and 3.4, namely, coj = and, u/y = (1 + xy) , hold. These statements 

cover the stability of the corresponding spectral and pseudospectral Jacobi approximations 
with variable coefficients. 

2. We should highlight the fact the stability assertion stated in Theorem 6.2 depends 
solely on the uniform bound of ct(xj) but otherwise is independent of the smoothness of a(x). 

3. The proof of Theorem 6.2 applies mutatis mutandis to the case of variable coefficients 
with a = a(x, t). If a(xj, t ) are C' 1 -functions in the time variable then (6.20) is replaced by 

< (1 + Const. A0K-,f m )|U-, < = (1 + Xj) ’ 

a{xj, t m ) 

and stability follows. 

4. We conclude by noting that the CFL condition (6.14b) depends on the quantity 

max i<i<^ rather than the minimal gridsize, , as in the constant coefficient case 

(compare (3.12)). This amplifies our introductory remarks in the beginning of this section, 
which claim that the 0(N~ 2 ) stability restriction is essentially due to the size of the Sturm- 
Liouville eigenvalues, A = 0(iV 2 ). Indeed, the other portion of the CFL condition 
requiring 

* „ “(Zi) 

At ■ 2 max - v — < rt 0 

i <i<N 1 - Xj ~ ' 

guarantees the resolution of waves entering through the inflow boundary x = 1. In the 
constant coefficients case this resolution requires time-steps At of size — . However, when 

the inflow boundary is almost characteristic, i.e., when a(l) ~ 0, then the CFL condition 
is essentially independent of Axmjn, for (6.21) boils down to At ■ 2a'{\) < tj 0 . In purely 
outflow cases such as (6.3), the time-step is independent of any resolution requirement at 
the boundaries and we are left with the CFL condition (6.5) stated in Theorem 6.1. 
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